Ionization Suppression of Diatomic Molecules in Intense Midinfrared Laser Field 
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Diatomic molecules (e.g., O2) in intense laser field exhibit a peculiar suppressed ionization be- 
havior compared to their companion atoms. Several physical models have been proposed to account 
for this suppression while no consensus has been achieved. In this letter, we aim to clarify the 
underlying mechanisms behind this molecular ionization suppression. Experimental data recorded 
at midinfrared laser wavelength and its comparison with that at near-infrared wavelength revealed 
a peculiar wavelength and intensity dependence of the suppressed ionization of O2 with respect to 
its companion atom of Xe, while N2 behaves like a structureless atom. It is found that the S-matrix 
theory calculation can reproduce well the experimental observations and unambiguously identifies 
the significant role of two-center interference effect in the ionization suppression of O2. 
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Ionization of molecules in intense laser pulses plays 
a central role in understanding strong field molecular 
physics. Most molecular strong field processes, such as 
above-threshold ionization (ATI), high- harmonic gener- 
ation (HHG), double ionization (DI), and Coulomb ex- 
plosion (CE) are derived directly from this fundamental 
process. Compared to atoms, for which the ionization 
and related phenomena have been well understood 
molecular ionization exhibits a large variety of peculiar 
behaviors due to its structural complexity and the extra 
nuclear degrees of freedom [2| . 

One of the most fundamental differences between 
molecule and its companion atom, i.e., the atom with a 
comparable ionization potential, is represented by their 
relatively distinct ionization probability. For example, 
when subject to a TkSapphire laser pulse at ~ 800 nm, a 
strong suppression has been observed in ionization prob- 
ability of diatomic molecule O2 compared to the rare 
gas atom Xe, while no suppression is seen in the di- 
atomic molecule N2 compared to its companion atom Ar 



[H, H[. Several theoretical models, including the KFR 
(or Keldysh- Faisal- Reiss) 0, ||, multielectron screening 
or MO-ADK 8] models, have been proposed to ad- 
dress this issue. The KFR model predicts that the in- 
terference between ionizing wave packets emitted from 
the two distinct nuclear centers can lead to ionization 
suppression for molecules (e.g., O2) with antisymmet- 
ric electronic ground states. The multielectron screen- 
ing model introduces a charge-screening correction to the 
tunneling theory. In the MO-ADK model, the difference 
between atomic and molecular ionization is attributed 
to different asymptotic behaviors of their ground state 
wave functions [8] . Though all the proposed models can, 
in principle, account for the experimental observations, 
mostly at near-infrared wavelength of 800 nm, no con- 
sensus on the underlying mechanism has been achieved 
so far. In addition, various ab initio methods have been 
applied to study this problem and qualitative agreement 
with the experimental results has been achieved [9l4lll|. 
however, no clear physical mechanism can be identified. 
Recently, alignment dependence of molecular ionization 
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has been measured and the comparison with theories 12 1 
shed more light on the underlying molecule specific ef- 
fects. 

Further understanding and clarification of the mech- 
anism behind the distinct molecular ionization requires, 
from experimental point of view, the extension of the 
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measurements into other wavelength range than the 
solely 800 nm. Indeed, recent experiment performed at a 
shorter wavelength of 400 nm [13| exhibited a similar be- 
havior of O2, and, however, a rather disparate behavior of 
N2 . The latter is found to have a higher ionization proba- 
bility compared to that of Ar for linear polarization while 
the difference vanishes for circular polarization. This has 
been explained by the resonance enhancement, a charac- 
teristic of multiphoton ionization (MPI) process, in N2. 
On the other side, this MPI resonance effect may also 
contribute significantly to the atomic or molecular ion- 
ization at 800 nm [l4j , frustrating an explicit comparison 
with the theory and a clear identification of the mecha- 
nism. 

In this letter, we perform a comparison study of in- 
tense field ionization between diatomic molecules (i.e., 
N2 and O2) and their companion atoms (i.e., Ar and Xe) 
at a midinfrared wavelength of 2000 nm. In contrast to 
previous studies carried out at shorter laser wavelengths, 
our experiment at long wavelength ensures that the ion- 
ization process falls deeply within the tunnel ionization 
(TI) regime [15|, [f6[ . The comparison of the data with the 
theory provides a crucial clue to the physical mechanism 
behind molecular ionization in TI regime, which is of spe- 
cial importance to the emerging held of ultrafast imaging 
of molecular structure and dynamics, wherein the ionized 
electrons from the molecules can be employed as a tool in 
imaging of molecular orbitals [l7l - fl9| and probing nuclear 
dynamics with attosecond resolution [20| . 

In our experiments, wavelength-tunable mid-infrared 
femtosecond laser pulses are generated by an optical 
parametric amplifier (OPA, TOPAS-C, Light Conver- 
sion, Inc.) pumped by a commercial TkSapphire laser 
system (Legend, Coherent, Inc.). This OPA system has 
been described in detail elsewhere [2l[ . The pulse energy 
from OPA is variable, before focused into the vacuum 
chamber, by means of an achromatic half- wave plate fol- 
lowed by a polarizer. A standard time-of-flight (TOF) 
mass spectrometer is used to register the ion signal. By 
means of a cryopump, the base pressure of the spec- 
trometer is maintained below I0 -8 mbar. At the end 
of the spectrometer, ions are detected with a microchan- 
nel plate as a function of flight time. The ion signal is 
further amplified, discriminated, and sent to a multihit 
time digitizer to generate TOF mass spectra. Depending 
on the laser intensity, the data point in the ion yield plots 
is obtained by averaging over 10 4 up to 6 x 10 6 laser shots 
at each intensity to ensure a sufficiently high statistical 
accuracy. 

In Figs. [T] (a) and (b), we present the measured ion 
yields of singly charged O2 versus Xe and N2 versus Ar, 
respectively, using a linearly polarized light at a center 
wavelength of 2000 nm. For comparison, data recorded 
at Ti: Sapphire laser wavelength of 800 nm are also shown 
in the inset. Very similar to previous studies at only 800 
nm 0], our data show that N2 and Ar have a parallel 
ionization probabilities also at 2000 nm. While for O2, 
a significant suppression of ionization yield compared to 
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Figure 1: (Color online). Experimental single ionization 
yields of diatomic molecules (i.e., N2 and O2) and their com- 
panion atoms (i.e., Ar and Xe) as a function of laser intensity 
at 2000 nm. The corresponding data recorded at 800 nm are 
also shown in the inset for comparison. 



Xe is found for both 800 and 2000 nm fields. 

In order to compare quantitatively the relative ion 
yields of the molecules with respect to their companion 
atoms and to provide more information for the benefit 
of understanding the underlying mechanism, we plot the 
ratio of ionization yield of O2 with respect to Xe in Fig. 
EJa), and N2 with respect to Ar in Fig. [2jb), respec- 
tively. It is found that the ratio of N^/Ar + keeps al- 
most constant, to be around 1, suggesting that the N2 
always behaves like a structureless atom, irrespective of 
the laser intensity and wavelength. On the other side, the 
ratio of O J /Xe + shows a strong dependence on the laser 
intensity. The ratio gradually increases with increasing 
intensity in both 800 nm and 2000 nm cases. More inter- 
estingly, the ratio shows a clear wavelength dependence. 
Two ratios almost coincide at about the lowest inten- 
sity used in the experiment but the ratio of 2000 nm ap- 
parently increases faster with intensity than that of 800 
nm, which indicates that, for the same intensities above 
5 ~ 6 x I0 13 W/cm 2 , the longer the wavelength, the less 
pronounced ionization suppression of O2 compared to Xe. 

This apparent dependence of the ratio of O^ /Xe + on 
the laser wavelength is inconsistent with the MO-ADK 
model prediction |8[ since ADK-type tunneling formula- 
tion is based on a quasi-static approximation which will 
give an ionization rate independent of the wavelength. 
On the other side, the multielectron screening model pre- 
dicts that the suppression of the O2 compared to Xe 
should be more pronounced with the increase of the laser 
wavelength [3, EH , which is also in disagreement with the 
experimental observations. In contrast, we will see below 
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Figure 2: (Color online). Experimental ratios of ionization 
yield between molecules and their companion atoms for both 
800 and 2000 nm, as a function of laser intensity. 

that this dependence of the ratio of 0^"/Xe + on both 
the wavelength and the intensity can be well understood 
within the S-matrix formulation and is closely related to 
the destructive interference of the two subwaves of the 
ionizing electron emerging from the two atomic centers 
of O2 with ground state of antibonding symmetry. 

In the S-matrix theory, single ionization rate for an 
atom or a molecule in a linearly polarized laser field with 
a vector potential A(i) = A$e z cos(cji) is given by [H, [22[ 

W =2vN e J2 k N(U P ~ Nuj) 2 f dk N J%(kn • a, ^) 

x|(0 fejv (r)|^(r))| 2 . (1) 

Here N e denotes the number of equivalent electrons and 
iVo is the minimum number of the photons needed to 
ionize the target, — y // 2(Nu! — l v — U p ) and k^v rep- 
resents the momentum of the emitted electron, with I p 
and Up the ionization potential and ponderomotive en- 
ergy, respectively. Jjv is a generalized Bessel function 
with a = Aoe z /u>. 4>k N (^) denotes the plane wave func- 
tion and ipi(r) is the ground state wave function of the 
atom or molecule. In our calculation, the wave func- 
tion of atomic ground state is approximated by the out- 
most single electron orbital, while for molecules we use 
the linear combination of atomic orbitals to simulate the 
molecular orbital (LCAO-MO) approximately 0, |24j . 
Moreover, according to ab initio calculation using time- 
dependent density- functional theory in Ref. [llj . the 
highest-occupied molecular orbital (HOMO) (1%) domi- 
nates in the ionization process of O2 molecule while both 
HOMO (3cr 9 ) and HOMO-1 (1tt„) play important role 
for N2. Therefore, only HOMO is considered for O2 and 
both HOMO and HOMO-1 are considered for N2 in our 
calculation. 

The calculated ratio of single ionization rates for 
molecules O2 and N2 with respect to their companion 



Figure 3: (Color online). Ratios of ionization yield between 
molecules and their companion atoms for both 800 and 2000 
nm, calculated from the S-matrix theory, as a function of laser 
intensity. 

atoms Xe and Ar as a function of laser intensity at 800 
nm and 2000 nm are shown in Figs. [3]Ja) and (b), re- 
spectively. It is found that the theoretical results are well 
consistent with the experimental data 25] . The ratios for 
/Ar + keeps around 1, irrespective of the field intensity 
and wavelength. In contrast, the ratios for O^ /Xe + show 
strong suppression and raise with the increase of the laser 
intensity. Moreover, the ratio of 2000 nm keeps higher 
than that of 800 nm for the whole intensity regime con- 
sidered. Especially, both experimental and theoretical 
data agree on that the ratio of O^ /Xe + depends linearly 
on the laser intensity (i.e., oc I), e.g., the ratio of 800 
nm increases by about 4 times when the laser intensity 
increases by about 4 times (see Fig. EJJa)). It is notewor- 
thy that this amount of increase is significantly larger 
than that predicted by the MO-ADK theory which gives 
a scaling of J 1 / 2 @]. 

Further comprehension of the distinct intensity and 
wavelength dependence of the ratio of 02VXe + can be 
drawn from the S-matrix formula, in which the single 
ionization rate of molecules consists of two major parts 
that comes from the Fourier transform of the ground 
state. For O2, the Fourier transform has the form 
^i(kjv) = C<I> at (kAr) sin(kiv -R/2) with C being the nor- 
malization factor. Here $ at (kjv) is the atom-like part, 
which is determined mainly by the atomic orbital consti- 
tuting the molecular orbital (2p x wavefunction for O2), 
and the trigonometric part sin(kjv • R/2) depends on the 
molecular structure and is associated with the interfer- 
ence effect between the wavepackets of the ionizing elec- 
trons centered at the individual nucleus Q . This trigono- 
metric part leads to the suppression effect since it always 
gives destructive interference when k^ ■ R <C it which 
is usually satisfied in the conditions considered in this 
work. 

This suppression effect can be clearly seen in Fig. 0Ja) 
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which shows the calculated ion yield ratio between the 
two cases with and without the trigonometric term in- 
cluded in Eq. ([I]), denoted by 0^" and O^*, respectively, 
for both 800 nm and 2000 nm laser wavelengths. When 
the trigonometric term is included, the ionization yield 
is significantly suppressed. Moreover, the ratio ascends 
with intensity and increases with laser wavelength. To 
understand these effects, we plot the momentum spectra 
of photoelectrons emitted from the ionization of molec- 
ular O2, calculated without the trigonometric term, for 
different laser intensities and wavelengths. As shown in 
Figs. 0Jb) and (c), calculated at 800 and 2000 nm, re- 
spectively, the contribution from large fcjv becomes more 
significant as the intensity increases. This will cause an 
increase of the interference term sin(kN • R/2) and as 
a consequence, the ascending ratio with the increase of 
the laser intensity, in agreement with the results in Fig. 
|4ja). In addition, for a certain laser intensity the mo- 
mentum spectrum becomes broader when the wavelength 
increases (see Fig. Scl)). Similar to the intensity effect, 
this results in an increasing ratio and less pronounced 
suppression effect for longer wavelength, which is consis- 
tent well with the experiment. Note that in Figs. HJb), 
(c) and (d) , for simplicity, we only show the electron dis- 
tribution for a specific molecular alignment angle of 45° 
with respect to field direction, since the ionization rate of 
O2 reaches maximum at about this angle [19|. It is wor- 
thy to mention that the distributions are not sensitive to 
the alignment angle in our calculations. Therefore, the 
interference effect from the trigonometric term originat- 
ing from the multi-center feature of molecules plays an es- 
sential role in the ionization suppression of O2 molecule. 
For N2 molecule, the situation is much more complicated. 
Besides the HOMO orbital, the HOMO-1 orbital con- 
tributes to the total ion yield and is becoming especially 
important at high intensities. Moreover, the HOMO or- 
bital is an admixture of both atomic s and p orbitals. 
The contributions of these two types of orbitals, which 
possess trigonometric terms of both cos(k/v ■ R/2) and 
sin(kAr • R/2), add coherently in the calculation [2g|. As 
a result, N^~/Ar + shows no suppression and is hardly de- 
pendent on the laser wavelength. 

Note that the ratio of Oj /Xe + has a scaling of oc I de- 
pendence on the laser intensity (see Figs. EJa) andlSJa)), 
in contrast to that the ratio shown in Fig. EJa), intro- 
duced solely by the interference effect, gives a scaling of 
oc I 1 / 2 . The other contribution accounting for the addi- 
tional I 1 / 2 comes from the effect of the atomic orbitals. 
It is noteworthy that the outmost orbital of O2 molecule 
is composed of a 2p x orbital while the outmost orbital of 
xenon atom is 5p z . Our calculation shows that the ratio 
between the ionization yield from these two orbitals will 
also increase by about 3^4 times when the intensity in- 
creases by about one order of magnitude, which gives an 
additional J 1 / 2 dependence of the ratio O^/Xe" 1 ". There- 
fore, the two-center interference effect, together with the 
atomic orbital effect, gives rise to the experimentally ob- 
served intensity dependence (oc I) of the 2 f /Xe + . 
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Figure 4: (Color online), (a): Calculated ratio between ion- 
ization yield of O2 with (denoted as Oj) and without the 
trigonometric term (denoted as O^"*) (see text), (b), (c) and 
(d): Normalized momentum spectrum of the photoelectron 
from O2 calculated without the trigonometric term at molec- 
ular alignment 45° with respect to the laser field. The field 
parameters used in the calculations are shown in the panels. 



A closer comparison between our theoretical simula- 
tion (Fig. EJa)) and experimental data (Fig. Ufa)) shows 
a perceptible discrepancy at very low intensity regime. 
The experimental ratio does not decrease with decreas- 
ing intensity while the theoretical ratio keeps dropping 
at about 4 ~ 5 x 10 13 W/cm 2 . This discrepancy may 
be partially attributed to resonance effect that becomes 
important when the laser intensity is low and the ion- 
ization is well in the multiphoton regime (l3| . In the 
multiphoton regime, the ionization channel via Freeman- 
resonance process [27| contributes significantly to the ion- 
ization yield. It is well-known that O2 molecule possesses 
more abundant highly excited states than Xe and hence 
provides more resonance channels for the ATI process, 
resulting in the increased ratio of 0^/Xe + comparing 
with that given by the S-matrix calculation in which all 
the resonance processes are ignored. 

In summary, we present a comparison study on the ion- 
ization of diatomic molecules (N2 and O2) and their com- 
panion atoms (Ar and Xe) at midinfrared wavelength. 
Our experimental data reveals that the ionization proba- 
bility of N2 is almost identical to its companion atom of 
Ar. In contrast, O2 exhibits a distinct suppression com- 
pared to Xe and more importantly, a strong dependence 
of this suppression on both the laser wavelength and in- 
tensity has been found. While this finding is in conflict 
with the molecular ADK formulation and multielectron 
screening model predictions, it can be well reproduced 



5 



by the S-matrix theory calculation, which considers the 
interference between ionizing wave packets emitted from 
the two ionic cores. This interference effect, together 
with the different intensity dependence of the ionization 
of atomic orbitals, accounts for the peculiar ionization 



behavior of O2 comparing to its companion atom Xe. 
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